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ABSTRACT 


This  is  a  signal  study  to  demonstrate  the  possibility 
of  reducing  near  surface  reverberations  due  to  geological 
effects  near  the  vertical  array  receivers.  The  signals  are 
deghosted  to  make  the  up-going  P-pulses  appear  similar  on 
all  of  the  vertical  array  sensors.  A  correlation  record  is 
computed  to  measure  the  similar  component  which  occurs 
jointly  on  all  of  the  array  elements.  The  coda  of  a  strong 
USC&GS  zero  focus  event  and  588  km  focus  event  were  consider¬ 
ably  simplified.  The  coda  of  a  USC&GS  60  km  event  from  the 
Aleutians  showed  sufficient  definition  to  improve  the  de¬ 
tection  of  pulses  occurring  after  the  first  P-pulse.  Ten 
unobservable  weak  signals  were  processed  and  three  were  de¬ 
tected,  based  on  proximity  to  Herrin  times,  amplitude  and 
character. 


1. 


INTRODUCTION 


The  concept  of  a  vertical  array  is  to  record  simultane¬ 
ously  on  several  transducers  which  are  stacked  in  a  deepwell 
which  is  3  km  or  more  in  depth.  We  are  interested  in  recording 
teleseismic  events  in  the  period  range  between  .5  to  1.5 
seconds.  This  is  equivalent  to  about  one-half  to  one-and-a- 
half  wavelengths  in  the  period  range  of  interest.  Ideally, 
we  want  to  place  instruments  at  one-fourth  wavelength  sepa¬ 
ration  at  the  shortest  wavelength  of  interest,  which  would  be 
approximately  .5  to  .7  km  apart.  For  our  data,  we  use  four 
deepwell  instruments  between  1*5  and  3  km  and  one  instrument 
near  the  surface.  This  is  a  good  configuration  based  on  the 
above  criteria.  The  instruments  are  run  on  the  same  gain 
Geotech  (1965),  which  makes  convenient  the  design  of  detectors 
based  on  the  combination  of  two  or  more  seismometers.  In 
addition,  new  developments  beyond  the  scope  of  this  report  in¬ 
clude  a  stacked  array  of  triaxial  seismometers  described  by 
Shappee  and  Douze  (1966). 

Although  there  is  little  background  literature  for  this 
kind  of  operation,  we  can  take  advantage  of  years  of  experience 
in  reflection  seismology,  where  vertical  stacking  of  an  ex¬ 
plosion  source  is  used  to  remove  distortion  of  near  surface 
layering,  mainly  the  ghost  reflection  from  the  surface.  By 
appi  priately  scaling  the  frequencies,  the  stacked  source 
array  in  reflection  prospecting  is  equivalent  to  a  maximum 
source  depth  of  .07  km  .  This  indicates  close  correspondence 
in  wavelengths  between  our  vertical  array  seismometer  experi- 


1 


ments  and  the  stacked  source  experiments  in  reflection  pros¬ 
pecting.  It  is  reasonable  to  anticipate  similar  results  and 
to  employ  similar  data  processing  schemes.  By  reciprocity, 
the  same  algorithms  developed  for  stacked  sources  can  be 
applied  to  stacked  receivers  by  reversing  the  sign  of  the  re¬ 
flection  coefficients.  Goupillaud  (1960)  and  Geotech  (1965) 
apply  inverse  filtering  to  remove  the  effect  >f  the  surface 
and  interbedded  reflections.  Various  authors  have  developed 
techniques  for  the  elimination  of  ghost  arrivals  interferring 
with  the  up-going  compressional  wave.  Among  these  are 
Schneider,  et  al.,  (1964),  Hammond  (1962)  and  Lindsay  (1960). 
The  processing  schemes  differ  in  the  amount  of  detail  re¬ 
quired  in  specifying  the  transmission  properties  over  the  span 
of  the  array,  in  the  required  control  of  source  intensity  or 
instrument  gain,  and  the  required  amount  of  time,  stability, 
and  accuracy  of  computer  processing.  To  design  a  feasible 
and  robust  processing  scheme  we  attempt  to  minimize  the  de¬ 
pendency  of  our  scheme  on  all  of  these  requirements.  The 
processing  is  designed  to  correct  strong  signal  measurements 
for  reverberation  effects  of  layers  near  the  surface  at  the 
receiver  depths.  We  also  investigated  the  possibility  of  im¬ 
proving  the  detection  of  weak  p-waves  through  processing 
vertical  array  measurements. 

1 . 1  Physical  Description  of  Array 

The  vertical  deepwell  array  at  Apache,  Oklahoma 
(Geographic  Coordinates  -  34°49'59.0"  N  and  98°26'09.0"  W) 
consists  of  five  Geotech  Model  11167  seismometers.  (The 
code  identification  and  operating  depth  of  each  instrument  is 


listed  alongside  of  its  respective  trace  in  figures  appearing 
this  report)  .  The  stratigraphic  and  velocity  pro¬ 
files  were  obtained  from  Geotech  (1964).  Basically,  the 
region  surrounding  the  AP-OK  array  consists  of  high-velocity 
limestones  (  v  «  6  km/sec)  overlying  an  igneous  basement  com¬ 
plex  having  velocities  in  the  neighborhood  of  5  km/sec. 

Field  measurements  on  selected  days,  e.g.,  during  June 
of  1965,  were  phased  for  an  up-going  P-pulse  by  recording  the 
data  with  magnetic  heads  adjusted  to  compensate  for  the  delay 
in  the  first  arrival  of  the  pulse.  Data  processing  was  de¬ 
signed  to  handle  both  kinds  of  recordings,  the  phased  and 
normal  un-phased  recordings. 

1 -2  Description  of  Data  Samples 

Seismic  data  was  selected  from  17  June  1965,  when 
instrument  phasing  was  done  and  from  15  July  1965,  in  which 
normal  non-phased  seismic  recordings  were  taken.  Deepwell  #5 
(DW5) ,  located  at  50  feet  below  the  earth's  surface,  was  used 
in  place  of  the  surface  trace.  It  resembled  very  closely  the 
waveform  at  the  surface  recorded  on  the  SPZ  instrument. 
Furthermore,  DW5  was  at  the  same  gain  settings  as  the  other 
deepwell  seismometers,  whereas  the  surface  trace  was  not. 

Samples  were  selected  to  demonstrate  events  having 
different  focal  depths  and  different  apparent  frequencies. 
Three  events  were  processed  with  all  aspects  of  the  processing 
displayed.  In  addition,  ten  other  small  events  were  processed 
to  demonstrate  any  large  increase  in  detection  capability  for 
a  vertical  array  processor.  For  these,  only  the  surface 
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measurement,  sum  trace,  and  correlation  traces  are  displayed 
for  comparison.  Table  1  briefly  describes  the  events  selec¬ 
ted  . 


1.3  Data  Processing 

Before  developing  the  deghosting  technique,  the 
deepwell  trace  was  simulated  from  a  surface  seismogram.  If 
a  deepwell  trace  can  be  simulated  by  adding  only  an  echo  from 
the  surface,  it  seems  entirely  plausible  that  by  removing  the 
ghost  reflection  from  the  deepwell  trace,  the  up-going  P- 
pulses  can  be  made  to  appear  similar  at  each  depth.  The  echo¬ 
time  delays  for  the  four  deepwell  instruments  are  illustrated 

in  Table  2.  These  were  computed  from  the  acoustic  log  verti- 

o 

cal  velocity  profile,  using  an  incidence  angle  of  15  ,  and 
were  adjusted  slightly  by  trial  and  error  to  obtain  the  best 
apparent  results. 


ISTRUMENT 

DEPTH  (Ft.) 

ECHO-TIME  DELAY 
( Sec  o ) 

DW1 

9454 

.95 

DW2 

7460 

o 

CO 

e 

DW3 

6465 

.,75 

DW4 

5456 

.55 

Table  2 
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Table  1.  Description  of  Events  Used  to  Measure  Vertical  Array  Teleseisraic  Signals 


Taking  K  as  the  number  of  digital  points  in  the 

echo,  assuming  a  trial  reflection  coefficient,  a,  (0.9  was 

used)  and  given  the  surface  or  near  surface  trace  X  ,  the 

i' 


phased  simulated  deepwell  trace  is 


Z.  -  X.  +  a  X. 
i  l  l  -  K 


The  normal  non-phased  simulated  deepwell  trace  is 
Zi  =  Xi  +  K/2  +  *  Xi  -  k/2  • 

For  simulation,  the  ghost  is  constructed  from  the  surface  re¬ 
flected  P-pulse .  The  surface  reflected  S-pulse  and  other 
reflected  pulses  are  neglected. 

1  *4  Segho sting  -  Method  I 

When  a  surface  or  near  surface  trace  is  available, 
the  deghosted  trace  can  be  constructed  by  simply  shifting  the 
surface  trace  by  K,  multiplying  by  ha,  and  subtracting  this 
resultant  from  either  the  deepwell  trace  (phased  vertical 
array)  or  the  deepwoll  trace  shifted  by  K/2  (non-phased  ver¬ 
tical  array).  For  phased  vertical  array  measurements,  the 
echo  or  ghost  is  removed  on  Z, 


2i  ’  Yi  -  *  X(i  -  K, 


and  for  non-phased  or  normal  deepwell  records 


Zi  Y(i  -  K/2 )  "  ^  X(i  -  K), 


-  5  - 


where  X,  Y,  and  Z  are  the  surface,  deepwell,  and  deghosted 
deepwell,  respectively. 

If  a  surface  trace  is  not  available,  the  ghost 
can  be  removed  by  an  inverse  filter,  as  developed  in  the 
next  section. 

1.5  Deghosting  Method  II 

It  would  seem  more  plausible  to  develop  a  de- 
ghosting  process  that  removes  the  echo  without  the  use  of  a 
surface  seismogram,  a  practical  technique  that  would  figur¬ 
atively  push  the  ghost  off  the  end  of  the  deepwell  trace. 

Given  the  reflection  coefficient  a,  and  the 

deepwell  trace  ^  the  phased  deghosted  trace  is  formed  as 

follows:  First  a  new  trace  Z,  ,  is  defined  to  be 

i,l 


Z,  _  Y t  _  —  X ,  .  _  \  •  OC. -  , 

1,1  1,1  (l  -  NT1)  1 

where  NT^  is  equal  to  the  product  of  the  sampling  rate  and  the 
echo  time  delay  (K)  .  This  above  equation  simply  shifts  and 

inverts  the  ghost  to  a  point  NT  farther  down  the  seismic 

1  2 
record  than  it  originally  was.  is  redefined  as  cl^  = 

and  NT  is  redefined  to  be,  NT^  =  2  .  NT^  and  a  new  trace 

Z.  .  is  formed  where,  Z,  _  =  X,  _  +  X, ,  .  .  a„  .  Here 

i,2  i,2  if2  (i-  NT2 )  2 

the  ghost  is  inverted  and  shifted  still  farther  (NT2  points) 
down  the  seismic  record.  This  iterative  process  is  continued 
until  the  ghost  reflection  is  pushed  off  the  seismic  record. 
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In  general 


Z. 

1 < 3 


“  X.  .  +  X 
1/1 


(i 


NT  .  ) 
3 


a . /  where  a  .  ,  .  =  a 
1  1  +  1  3 


a  , 
3 


NTj  +  ^  =  2  .  NT j ,  and  NT^  <  N  (The  number  of  data  points 


in  the  deepwell  trace).  (Refer  to  Appendix  I  for  complete 
discussion  of  Z-transform) . 


For  non-phased  deghosting,  the  same  technique,  as 
that  developed  above,  applies  except  the  final  deghosted  trace 
is  time-shifted  by  K/2  to  the  right  in  order  to  align  the 
deepwell  with  the  surface  trace. 

1 . 6  Vertical  Velocity  Filtering 

(a)  Phased  sum  of  vertical  array  -  The  phased  sum 

of  the  vertical  array  is  formed  by  merely  adding  all  traces 

point  by  point.  That  is 
N 


i=l 


where  i  is  an  index  which  numbers  the  sensors  and  N  is  the 

number  of  seismometers  in  the  set  of  vertical  arrays  X.  .  „ 

1  •  3 

(b)  Phased  sum  of  deghosted  vertical  array  - 
The  same  formula  applies  for  the  deghosted  phased  sum  ex¬ 
cept  X.  .is  the  set  of  deghosted  vertical  arrays. 

(c)  Multi-channel  minimum  variance  estimate  of 
a  signal  fixed  on  all  elements  of  the  vertical  array  -  The 
formula  for  the  correlation  trace  is, 
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X. 

1 


J 


=  < 


<  <xlX2>  x3> 


n 

j  =2 


<X,  .X. 

1-]  1 


.> 

3 


where  J  is  the  number  of  deepwell  channels  on  the  seismogram 
(refer  to  Appendix  2  for  derivation),  i  is  the  number  ordering 
the  sequence  of  points  on  the  time  axis,  and  the  symbol  <  > 
represents  smoothing  with  a  moving  taper  or  time  window. 


1.7  Discussion  of  Programs 


The  following  eight  digital  programs  are  available 
to  either  remove  the  ghost  reflection  or  simulate  a  deepwell 
trace: 

1*  PROGRAM  SSIMDEEP  works  on  a  surface  trace  to  simulate 
a  phased  deepwell  seismic  record.  Given  the  echo 
times  between  the  deepwell  and  the  surface  and  further 
given  the  reflection  coefficient,  the  phased  deepwell 
trace  is  computed  by  SUBROUTINE  SSIMGO.  The  sum 
traces  of  both  the  deep  and  simulated  deepwell  traces 
are  computed.  Plots  are  obtained  of  the  deepwell 
seismogram,  simulated  deepwell  seismogram,  and  the  two 
above  sum  traces. 


2*  PROGRAM  SGOSTSEI  removes  the  ghost  reflection  when  a 
surface  or  near  surface  trace  is  available  at  the 
same  gain.  Given  the  echo  times  and  the  reflection 
coefficient,  the  phased  deghosted  trace  is  computed 
by  SUBROUTINE  SDEGOSTE .  The  sum  and  correlation 
traces  are  also  computed.  Plots  are  obtained  of  each 
deepwell  and  simulated  deepwell  seismogram  and  their 
corresponding  sum  and  correlation  traces. 
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3.  PROGRAM  SDEGSEI  removes  ghost  reflection  on  a  deepwell 
seismogram  with  the  use  of  SUBROUTINE  SDEGOSTI,  given 
only  the  echo-time  and  the  reflection  coefficient. 

Plots  are  obtained  of  deepwell  and  deghosted  deepwell 
seismograms,  sum  traces  arid  correlation  traces. 

4-6.  PROGRAMS  SIMDEEP ,  GOSTSEIS,  and  DEGSEIS  work  in  place 
of  1,  2,  and  3,  respectively,  when  normal  non-phased 
seismic  data  is  used.  Calls  are  made  to  SUBROUTINES 
SIMGO ,  DEGOSTE,  and  DEG0STE1 . 

7-8.  PROGRAMS  SVARRAY  and  VARRAY  are  used  on  phased  and 
non-phased  data,  respectively,  in  place  of  SDEGSEI 
and  DEGSEIS  when  only  plots  of  surface,  deghosted 
sum,  and  correlation  traces  are  wanted.  A  bandpass 
filter  is  applied  to  the  input  traces. 

2 .  RESULTS 

The  events  were  selected  to  demonstrate  measurements 
of  signals  of  various  types.  The  results  of  array  processing 
are  shown  in  the  following  figures; 

I  shows  the  raw  data  measurements,  phased  sum,  and 
phased  correlation  traces 

II  compares  the  deepwell  measurements  to  simulated 
deepwell  measurement 
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III  shows  the  phased  deepwell  measurements  deghosted 
by  subtraction  of  the  surface  measurement  and 
includes  the  sum  and  correlation  traces 

IV  shows  the  phased  deepwell  measurements  deghosted 
by  means  of  an  inverse  operator  which  removes 
the  surface  echo  and  includes  the  sum  and  corre¬ 
lation  trace 

The  figure  numbers  are  followed  by  a  letter  (such  as  II-B); 
the  letter  designates  the  event*  The  events  ares 

A  an  Aleutian  event  with  an  extremely  complex  coda 
and  focus  of  about  60  km,  selected  to  demonstrate 
possible  removal  of  the  part  of  the  complexity 
due  to  station  reverberations  to  facilitate 

picking  of  multiple  transmissions  of  P  and  depth 
phases 

B  an  event  from  the  Philippines  with  a  focus  of 

about  600  km,  selected  to  demonstrate  improvement 
in  the  simplicity  and  reduction  of  false  multi¬ 
ples  in  the  coda 

C  a  simple  event  from  Kazakh,  with  USC&GS  focus 
at  the  surface 

In  figures  IA-IC,  the  apparent  contrast  between  signal 
and  noise  appears  to  be  greater  at  the  surface  than  in  the 
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Figure  I-A.  Vertical  Component  Vertical  Array  Measurements— Aleutians 
USCGS  Focus  67  km,  Magnitude  5.2,  Date  6/17/65,  Time  19:05:9.1 
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Philippine  Islands,  USCGS  Focus  588  km. 
Magnitude  5.8,  Date  7/15/65.  Time  18:33:29.9 
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deepwell;  the  6  second  and  .5  second  noise  is  accentuated  in 
the  deepwell  measurements,,  Based  on  these  observations,  AP-OK 
does  not  appear  to  be  the  best  possible  selection  to  show  im¬ 
provements  by  vertical  array  processing.  All  of  the  sum  traces, 
phased  for  an  upgoing  vertical  P-wave,  appear  to  yield  higher 
ambient  and  signal  generated  noise  in  contrast  with  the  ap¬ 
parently  visible  signal  phases.  The  correlation  trace  shows 
moderate  improvemer. u.  in  signal  contrast,  with  better  definition 
of  multiple  arrivals  and  depth  phases  in  Figure  I  A  and  con¬ 
siderable  cleanup  of  the  coda  in  Figures  I  B  and  I  C„ 

Figures  II  A  -  II  C  demonstrate  that  the  deepwell 
seismogram  can  be  simulated  from  the  surface  trace  by  only 
considering  the  free  surface  reflection,  and  neglecting  inter- 
bedded  reflections.  Since,  for  the  forward  problem  of  de¬ 
riving  the  deepwell  seismogram  from  the  surface  seismogram, 
it  is  adequate  to  consider  only  the  P  surface  reflection,  it 
then  follows  that  in  designing  the  inverse  operator  it  is  also 
adequate  to  consider  only  the  free  surface  reflection.  We 
note  that  the  correspondence  between  deepwell  measurements  and 
the  simulations  is  generally  quite  good,  although  a  smaller 
uncorrelated  component  can  be  observed  which  can  be  attributed 
to  interbedded  reflections,  scattering,  etc.  This  component 
should  also  be  uncorrelated  between  deepwell  sensors  and  at 
least  partially  removed  by  the  correlation  trace. 

Figures  III  A  -  III  C  and  IV  A  -  IV  C  allow  us  to  compare 
two  different  methods  for  removing  the  surface  echo.  Method  I 
shown  on  Figure  III  is  conceptually  simple  and  involves  only 
subtracting  the  shifted  surface  trace  to  remove  the  echo. 
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Figure  II-B.  Comparison  of  Deepwell  and  Simulated  Deepwell  seismoc 
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Figure  II-C.  Comparison  of  Deepwell  and  Simulated  Deepwell  Seismograms 
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Figure  IV-B.  Deghoated  Vertical  Array  Measurements; 
Method  II  (Reflection  Coefficient  0.7  on  all  traces) 
Earthquake  Philippine  Islands,  USCGS  Focus  588  km, 
Magnitude  5.8,  Date  7/15/65,  Time  18:33:29.9 


Method  II  shown  on  Figure  IV  utilizes  an  inverse  operator  to 
remove  the  echo.  Both  methods  do  an  adequate  job;  the  corre¬ 
lation  trace  for  Method  II  has  a  slightly  cleaner  coda  and 
this  method  has  the  advantage  of  not  requiring  a  quiet 
surface  trace  and  exact  control  of  instrument  gain.  Since  the 
computation  time  of  Method  II  is  extremely  fast,  it  is  practi¬ 
cal  to  utilize  it  for  production  computations  of  the  corre¬ 
lation  trace. 

Figures  V,  VI,  and  VII  show  a  sample  of  ten  events  for 
which  we  seek  a  preliminary  demonstration  of  the  detection 
capability  of  this  vertical  array.  The  events  are  unobservable 
with  Herrin  times  indicated  by  an  arrow.  Figure  V  shows  the 
surface  measurements,  Figure  Vi  the  deepwell  deghosted  phased 
sum  measurements,  and  Figure  VII  the  deghosted  correlation 
trace.  All  of  the  abo\ e  outputs  are  filtered  with  a  broadband 
filter  set  at  a  center  frequency  of  1.25  cps .  We  consider  the 
detection  a  success  only  if  a  phase  appears  clearly  visible. 

The  fraction  of  total  events  detected  is  10%  for  the  filtered 
surface  trace,  0%  for  the  phased  sum,  and  30%  for  the  corre¬ 
lation  trace. 

3.  CONCLUSIONS 

Vertical  array  processing,  based  on  previous  experience 
with  stacked  sources,  appears,  as  expected,  to  be  an  effective 
means  of  reducing  the  reverberations  and  complexity  of  the 
coda  caused  by  dear  surface  reverberations  at  the  receiver. 

As  to  increasing  the  detectability  of  weak  P-pulses  based  on 
nearly  isotropic  vertical  velocity  filtering,  the  preliminary 
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Figure  V-l.  Surface  Vertical  Components  (Filtered) 
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results  suggest  some  improvement  is  probable .  Both  deghosting 
methods  work  well,  the  inverse  deghost  mg  operator  favored 
as  it  does  not  require  exact  control  of  the  gam  on  all  instru¬ 
ments  and  does  not  require  a  quiet  surface  recording.  The 
correlation  traces  used  for  detecting  the  upgoing  P-pulses 
impose  the  strongest  possible  requirement  that  the  signal  be 
fixed  jointly  on  all  channels.  For  zero  lags,  the  processing 
time  is  approximately  3/4  minutes  per  minute  of  data  for  all 
five  elements  of  the  vertical  array.  Similarly,  using  lags, 

the  processing  time  is  approximately  lh  minutes  per  minute  of 
data . 
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appendix  I 


Deghostincr,  Methnd  tt 

First,  obtain  the  forward  operation  2-transform  which 
generates  the  deepwell  trace  from  the  surface  trace 

D(t>  =  *s  sitt^lusit-^i 

The  reflection  coefficient  i , 

v.f  and  two-way  echo  time  is  t 

The  basic  shift-time  unit  for  purpose  of  ^  ° 

t- -<rpose  of  the  Z-transform  is 


t  . 
o 


-k 


D(Z)  -  h  (z  2  +  a  z+is)  s  (2) 


h  s(z)  =  — — 

1  +  m  z 


‘  [z*2  D  (2)  ]  [l  -  a  Z  +  a2  z2  .  . J 


H 


The  2 2  means  shift  one  half  the  echo  time  before  or  after 
performing  the  right  hand  infinite  sequence  of  operations. 
These  operations  involve  successive  scale  factoring 
shifting,  and  adding  or  subtracting  as  indicated.  The  sta¬ 
bility  of  the  operation  under  noise  and  speed  of  convergence 
is  improved  by  re-arrangement  of  the  operations  as  follows: 

[l-*2+*2  22  -....]=[(!  -a2,  +a2z2(1_az) 

+  h4  24  {(i  -  a  Z)  +  a2  z2  (1  -  a  2) }  +  ....] 
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Note  that  each  iterative  shift  operation  on  the  right  side 
doubles  the  interval  between  the  ghost  reflection  and  primary 
rather  than  increasing  it  linearly  as  on  the  left  side  so 
t.'nat ,  term  by  term,  the  recursive  formula  on  the  right  side 
converges  much  faster.  The  factor  of  gain  in  processing 
time  is  (log2N)/N  or  about  100  for  a  one-minute  record  at 
20  points/second.  In  the  time  domain,  the  above  operator  can 
be  expressed  as  a  recursive  filter. 

0  (t)  =  D  ( t )  -  nr.  D  (t  -  tQ) 

0  +1(t)  =  0.(t)  +  nrkD(t  -  ktQ) 

k  =  2 ^  (j  ~  1,2, . . . . »J ) 


and  J  t  <  T 
o 


where  T  is  total  record  length  in  seconds. 


>5  S(t  +  to/2)  -  0j(t) 
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APPENDIX  II 


Estimate  of  the  Jointly  Correlated  Signal  from  Multi-Channel 


The  minimum  variance  estimate  of  the  co  -linear  component 
on  two  time  series  for  complex  signals  is  taken  from  Mason  and 
Zimmerman  (1960).  The  energy  is  defined  as  the  squared  magni¬ 
tude  of  the  signal.  The  correlation  coefficient  is'  a  complex 
number  designating  the  amplitude  ratio  and  phase  (or  time) 
shirt  between  the  two  signals.  For  the  integrated  square  of 
the  residuals,  we  find  the  energy  integral 


Re  designates  real  part  and  asterisk,  complex  conjugation. 
The  middle  integral  on  the  right  is  a  complex  number  and  may 
be  expressed  in  polar  form 


J 


V. 


dt 


A  l 


j9 


We  adjust  the  magnitude  and  polar  angle  of  c12  to  minimize  the 

energy  integral.  The  optimum  angle  for  C  is  0 

12 


Substituting  in  the  energy  integral 


2 


A 
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We  minimize  the 


energy  by  adjustment  of  the  magnitude  of  C 


!—  J|V1  -  C12  V2\  2  dt  -  -  2  A  ♦  2|  Cu|  J|v2|2  dt  -  0 


2  dt 

0 

Multiplying  both  sides  by 

r 

J  v1  V  dt 

°12  ’  n  12 

J|v2|  dt 

If  the  signals  are  pulses  imbedded  in  noise  on  the  two  channels 
of  pulsewidth  T,  we  may  choose  to  replace  the  above  integration 
by  an  average  or  summation  restricted  to  a  time  window  T, 
symbolized  as  <  >  .  For  this  case  we  use  a  moving  time-window 

for  estimating  C^2  at  each  point  along  the  time  series 


C 


12 


(t) 


+T/2 

J  v^t  '  -  t)  v2  (t  '  -  t)  dt  7 

-T/2 _ 

+T/2 

J  |v2(t'  -  t)|  2  dt' 

-T/2 


For  a  signal  which  is  aligned  or  linearly  polarized  on  the  two 

channels,  we  can  take  the  signal  as  real  in  which  case 

v  =  v  For  case  where  the  signal  is  band-limited  and 
2  2 
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imbedded  in  broad-band  noise  we  compute  the  correlation  coef 
ficient  with  lags. 


If  the  signal  is  both  band-limited  and  linearly  polarized  on 
the  two  components, 


+T/2 


Generally,  a  better  estimate  of  an  imbedded  signal  can  be  ob¬ 
tained  if  the  co-ordinates  v  and  v  are  rotated  to  v  '  and 

12  1 

V2  .  The  rotation  which  maximizes  S/N  ratio  is  given  by 
Sax  (1966).  The  operator  C12  is  used  as  a  time  varying  filter 
to  extract  an  estimate  of  v.^  from  v2 
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or  it  can  be  used  as  a  gain  control  voltage  to  estimate  the 
imbedded  signal. 


A 

V 


l(t) 


‘  Cl2(t)  V2(t) 


If  preferred,  the  Wiener  criteria  is  used  to  estimate  the 
signal 


vi(t>  ' 


*  h  V2  > 


(<  V. 


<  V. 


2  >  >* 


v-^t) 


and  similarly  as  before  with  lags,  we  replace  products  in  the 
numerator  by  convolutions. 

For  the  multi-channel  case,  we  have  a  fixed  signal  im¬ 
bedded  in  N  noise  channels.  We  ask  that  the  signal  be  jointly 
correlated  on  all  channels.  As  a  fast  method  for  measuring 
the  jointly  correlated  signal,  we  make  a  chained  estimation  of 
the  v1(t)  on  the  N  channels  as  follows  for  the  zero  lag  case. 

The  minimum  variance  estimate  of  from  is 


1,2 


<  V1  V2  > 

(t)  *  - -  v- 


of  v  from  v  and  v_  is 
1  2  3 


ft 


Vl, (2,3) 


A 

<V 


1,2 


V  V3 


<  <vt  V  >  V  > 
1  2  3 


<  vj  >  <  V 


V 


3 


and  v1  from  ,  v  , 


N 


1,(2, 


,  N) 


<  •••  <  <V  V  >  V  „>  ...  >  V  > 

■ - 1 _ 2_  3  N 


A  1 2  | 

<  v„  >  <  V 

2 

<  I2  V* 

1  2I  | 

3 

'*  <  vn|  > 

and  with  lags  (including  zero  lag  as  a  special  case) 


<  • 

•  •  t 

e  wi 

®v  ®vv 

>  .  . 

••••  ®  V 

< 

V2 

V 

A 

CN 

V3 

2 

...  < 

1 

VN 

j2>  VN 

This  is  one  of  the  operators  used  in  this  report.  For  the 
case  of  lags  we  replace  the  products  in  the  numerator  by  con¬ 
volutions,  but  practical  time  of  computation  considerations 
alone  rule  out  anything  but  a  small  number  of  lags.  Note  in 
the  above  formula  that  v,  is  estimated  from  vn  which  may  be 
selected  from  any  of  the  available  channels.  The  correlation 
is  most  useful  .  f  all  of  the  channels  are  approximately  the 
same  quality,  and  the  channel  selected  for  vN  should  be  the 
best  available;  in  our  case,  this  is  the  surface  seismogram. 
If  the  seismograms  are  of  widely  varying  quality,  the  above 
formula  should  be  generalized  to  include  rotation  and  a  test 
to  eliminate  bad  channels.  Since  by  hypothesis  we  require 
correlation  on  all  channels,  we  also  observe  that  it  is  not 
necessary  to  consider  correlations  on  all  possible  pairs  of 
sensors,  resulting  in  considerable  time  savings  in  the  calcu- 


lation  of  the  correlation  trace  The  correlation  (2  ,N) 

may  be  generalized  to  include  lags. 


lation  of  the  correlation  trace.  The  correlation  (2 , „ „ „ „ ,N) 
may  be  generalized  to  include  lags. 
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It  ABSTRACT 


This  is  a  signal  study  to  demonstrate  the  possibility  of 
reducing  near  surface  reverberations  du<2  to  geological  effects 
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first  P-pulse.  Ten  unobservable  weak  signals  were  processed  and 
three  others  were  detected,  based  on  proximity  to  Herrin  times, 
amplitude  and  character. 

■■iMMMMBHHHMMMMMMMMHHHr  .VBH 

DD  1473  Unclassified 

Security  Classification 


Unclassified 


Security  Classification 


LINK  A 

LINK  e 

LINK  C 

1  ROLE 

•T 

ROLE 

„T 

ROLE 

WT 

Vertical  Seismic  Array  Recording 
Deepwell  Studies 
Seismic  Instrumentation 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY;  Enter  the  name  end  addreaa 
of  the  contractor,  tube ontrac tor,  grantee,  Department  of  De- 
fenae  activity  or  other  organization  ('corporate  author)  issuing 
the  report. 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over¬ 
all  security  claasification  of  the  report.  Indicate  whether 
“Restricted  Data*'  ia  included.  Marking  is  to  be  in  accord- 
1  V.CS  witN  expropriate  security  regulations. 

26.  GROUP:  Automatic  downgrading  la  specified  In  DoD  Di- 
r'~*rtive  5200. 10  and  Armed  Forces  Industrial  Manual,  Enter 
|  ...  group  number,  Alao,  when  applicable,  show  that  optional 
l.kings  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ed. 

3.  REPORT  TITLE:  Enter  Lite  complete  report  title  in  all 
capital  lettera.  Titles  in  all  cases  should  be  unclassified. 

If  s  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  in  ail  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g. ,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  ia 
covered. 

5.  AUTHOR(S):  Enter  the  name(s)  of  authors)  as  shown  on 
or  in  the  report.  Entei  laat  name,  first  name,  middle  initial. 

If  military,  nhow  rank  and  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report  as  day, 
month,  yea/;  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7  a.  TOTAL  NUMBER  CF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  I.e,,  enter  the 
number  of  p«^;**a  containing  information. 

7b.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  numher  of  the  contract  or  grant  under  which 
he  report  was  written. 

116,  8c,  U  8 d.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  aa  project  number, 
Subproject  number,  system  numbers,  tsak  number,  etc. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  t?  e  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  muat 
be  uni  qua  to  this  report. 

96,  OTHER  REPORT  NUMBER($):  If  the  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  sponsor),  alao  enter  this  number(s). 

10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those. 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  “Qualified  requestera  may  obtain  copies  of  this 
report  from  DDC" 

(2)  "Foreign  announcement  and  diaaemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
uaera  shall  request  through 


(4)  “U.  S.  military  agencies  may  obtain  copies  of  tKa 

report  directly  from  DDC  Other  qualified  uaerr 
shall  request  through 


(5)  "All  distribution  of  this  report  ia  controlled.  Qual¬ 
ified  DDC  uaera  shall  request  through 


If  the  report  haa  been  furnished  tc  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  aale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

11.  SUPPLEMENTARY  NOTES:  Uae  fnr  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  aponaoring  (pay¬ 
ing  for)  the  research  and  development.  Include  addreas. 

13-  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
aummary  of  the  document  indicative  of  the  report,  even  though 
it  may  alao  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port.  If  additional  space  is  required,  a  contin'  tion  sheet  shall' 
be  attached. 

It  is  highly  deairable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in¬ 
formation  i.i  the  paragraph,  represented  aa  (TS).  ( S ).  (C),  or  (U). 

There  ia  no  limitation  cn  the  length  of  the  abatract.  How¬ 
ever,  the  suggested  lengtl'  ia  from  150  to  225  worda. 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  snd  may  be  used  aa 
index  entries  for  cataloging  the  report.  Key  words  muat  be 
selected  »o  that  no  security  claasification  la  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  sn  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional. 


GTO  886*551 


Unclassified 

Security  Classification 


